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ABSTRACT. The pivotal role of electrical energy in propelling a nation’s economic development can- 
not be overstated. A deficiency in the efficient consumption of electricity not only hinders economic 
growth but also constitutes a potential threat to national security. Recognizing the gravity of this rela- 
tionship, scholars and policymakers have directed substantial attention towards exploring electricity 
consumption patterns and developing predictive models. Consequently, many models and methodolo- 
gies have been employed to forecast electrical energy consumption. 

Within this context, the present research delves into a comparative analysis of two prominent fore- 
casting methods, ARIMA and GM (1,1) grey modelling techniques, specifically focusing on predicting 
electricity consumption in Algeria. The assessment of these models’ performances was conducted us- 
ing the average percentage of absolute MAPE, employing annual data collected from Algeria spanning 
the extensive period from 1982 to 2020. 

The research findings underscore the paramount importance of accurate forecasting in this domain. 
Notably, the ARIMA model (1,1,0) emerged as the frontrunner, exhibiting superior predictive capabilities 
when juxtaposed with the GM (1,1) model, as evidenced by the MAPE standard. This nuanced examina- 
tion contributes to the scholarly discourse on electricity consumption prediction, offering insights that 
can inform strategic decision-making and policy formulation in pursuing sustainable and secure energy 
practices. 


KEYWORDS: FORECASTING, ELECTRICITY CONSUMPTION, AUTOREGRESSIVE INTEGRATED MOVING 
AVERAGE (ARIMA) MODEL, GREY MODELLING (1,1). 
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1. INTRODUCTION 


Electrical energy is one of the main pillars in 
the economic development of countries, as any 
lack of consumption of this energy may lead to 
a slowdown in economic development and eco- 
nomic growth and may even threaten national 
security (Panklib et al., 2015). [1] Therefore, elec- 
tricity consumption has become a topic of great 
importance because of the increasing demand 
for energy worldwide, mainly fueled by grow- 
ing economic activities (Kandananond, 2011). [2] 
Given this fact, a sound forecasting technique 
is essential for accurate investment planning of 
energy production, generation, and distribution 
(Bianco et al., 2010). [3] 

In recent decades, many methods in the liter- 
ature for forecasting used different techniques for 
electricity consumption prediction. The most com- 
mon models for this purpose are the autoregres- 
sive integrated moving average, artificial neural 
networks (ANN), support vector regression (SVR), 
hybrid models, and fuzzy logic methods. 

Algeria has recently worked to develop the 
state in all fields, whether economic, social, po- 
litical, etc. 

As the Electric energy index indicates a factor 
for the development of any country, the govern- 
ment of Algeria should monitor electricity con- 
sumption growth; forecasting is essential for the 
practical application of energy. Accurate forecasts 
of electricity consumption are vital when demand 
is growing faster. On the other hand, the values of 
electricity consumption in Algeria can be viewed 
as volatile and increasing. 

The present paper aims to analyze and fore- 
cast electricity consumption in Algeria utilizing 
The grey forecasting GM (1,1) model, invented by 
Deng, 1982 (Wang et al., 2018) [4], which has been 
regarded as one of most widely utilized prediction 
approaches because its requirement of few sam- 
ples and high forecasting accuracy, and autore- 
gressive integrated moving average(ARIMA). 

The rest of the paper is organized as follows. 
In Section 2, some related works are introduced. 
Section 3 describes the methodology. Section 4 
presents the main results and discussion. Finally, 
the conclusions and future works are presented in 
Section 5. 
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2. LITERATURE REVIEW 


Electricity consumption forecasting has re- 
ceived prominent attention among researchers 
today using different models (Ogcu et al., 2012). 
[5] This study employed seasonal SVR and ANN 
models to develop the best model for predict- 
ing electricity output between 2010 and 2011. The 
study’s results provided empirical evidence that 
electricity consumption in Turkey will increase 
soon. (Rathnayaka, et al 2014). [6] GM (1, 1) and 
Grey Model (1, 1) have been used for forecasting 
results. The performance of the proposed tech- 
nique has been compared with the existing Au- 
to-regressive moving average forecasting model. 
Annual power generation and forecasting data in 
Sri Lanka were used. MAPE, MSD and MSE accu- 
racy testing results show that GM (1, 1) outper- 
formed compared to model fitting and forecast- 
ing. (Hussain et al, 2016). [7] This paper applies 
Holt-Winter and Autoregressive Integrated Mov- 
ing Average (ARIMA) models on time series sec- 
ondary data from 1980 to 2011 to forecast total 
and component-wise electricity consumption 
in Pakistan. Results reveal that Holt-Winter is 
the appropriate model for forecasting electric- 
ity consumption in Pakistan. (SISMAN, 2017), [8] 
proposition of a new approach by comparing 
grey prediction and ARIMA models with the Mod- 
el of Analysis of the Energy Demand in Turkey 
from 1970 until 2013. Results have revealed that 
ARIMA and GP give better results than MAED to 
forecast long-term perspectives. In addition, it 
appears that the original ARIMA and GM (1,1) mod- 
els have powerful forecasting models. (Amber et 
al., 2018). [9] This paper uses Multiple Regression 
(MR), Genetic Programming (GP), Artificial Neural 
Networks (ANN), Deep Neural Networks (DNN) 
and Support Vector Machine (SVM) to forecast 
the electricity consumption of an administra- 
tion building located in London, United Kingdom. 
Results demonstrate that ANN performs better 
than all other four techniques. (Jain et al., 2018). 
[10] In this paper, an attempt is made to forecast 
electricity consumption using the ARIMA mod- 
el. The mean absolute percentage error (MAPE) 
measures forecast accuracy. Results show that 
the ARIMA model has the potential to compete 
with existing techniques for electricity consump- 
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tion forecast. (Leite et al., 2022). [11] This paper 
aimed to apply time-series forecasting models 
(the Holt-Winters, SARIMA, Dynamic Linear Mod- 
el, and TBATS (Trigonometric Box-Cox transform, 
ARMA errors, Trend, and Seasonal components 
and NNAR (neural network autoregression) and 
MLP (multilayer perceptron). The results indicate 
that the MLP model obtained the best forecast- 
ing performance for the electricity consumption 
of the Brazilian industry under analysis. (Bilgili 
et al, 2022). [12] In this study, four time-series 
methods: long short-term memory (LSTM) neu- 
ral network, adaptive neuro-fuzzy inference sys- 
tem (ANFIS) with subtractive clustering (SC), AN- 
FIS with fuzzy means (FCM), and ANFIS with grid 
partition (GP) were implemented for forecasting 
electricity consumption. Results reveal the LSTM 
model generally outperformed all ANFIS models. 


3. METHODOLOGY 
3.1. The ARIMA Model 


The model was proposed by Box-Jenkins in 
1970, as It consists of autoregressive models, which 
are denoted by (AR) of order (p), and moving av- 
erage models, which are denoted by (MA) of order 
(q), and (d) indicates the number of times the time 
series differences should be taken to make it sta- 
ble (In most cases d = 1) (Fattah et al., 2018). [13] 

1- TheAutoregressive (AR) Model 

The equation below takes the following: 
(Goshit, &lorember, 2018) [14] 


Ve = A+ OyyYp_1 + O2Vt-2 + °° + OyVe_-p + & 


2- TheMovingAverage (MA) Model 
can also be designed as follows: 


Ve = UF Et + PxEe_1 + P2Et-2 + + Dg Et—q 


3- TheAutoregressivelntegratedMovingAver- 
age (ARIMA) 
The formula for ARIMA could be represented 
as: (Grigonyte&Butkevicilité, 2016) [15] 


Ve = A+ OVe-1 + O2Vt-2 + °° + OnVt-p + 
+ Ep + PyEe_1 + P2Ep-2 + °° + PgEt—q 


GLOBALIZATION AND BUSINESS #16, 2023 


Developing the ARIMA model requires to fol- 
low four basic steps: identification, estimation, 
diagnostic checking, and forecasting (Bari et al., 
2015) [16] 

differencing to achieve stationary Stability 
analysis of the time series is necessary for the 
modelling process. We consider a series stable if 
its arithmetic average, variance, and covariance 
over time are constant (ATIL& MAHFOUD, 2021). 
[17] We can find out, however, that the series is 
unstable by The autocorrelation functions (ACF) or 
partial autocorrelation functions (PACF) (Udom&- 
Phumchusri, 2014) [18], or Augmented Dicky-Fuller 
Test’s (ADF) and PP Test’s. 

Step one is to identify (p, q), which the orders 
of the AR model and MA model and is the most im- 
portant, using the autocorrelation (ACF) and par- 
tial autocorrelation (PACF) functions are taken to 
identify the models because they show different 
features of the functions (Guha&Bandyopadhyay 
2016). [19] For AR (p), the ACF tails off at the or- 
der of p but PACF cutoff; for MA (q), the ACF cut- 
off but PACF tails off at the order of q; for ARMA 
(p, q), none of the ACF and PACF tail off, and to 
select as the best suitable model for forecasting 
(Mondal et al, 2014) [20], using BIC (Bayesian In- 
formation Criterion) and AIC (Akaike Information 
Criterion) values. Step two is Estimating the co- 
efficients of the formulation using the Maximum 
Likelihood method (Eissa, 2020). [21] Step three is 
the Checking of the Model, where the model must 
be checked for adequacy by considering the prop- 
erties of the residuals and whether the residuals 
from an ARIMA model must have a normal distri- 
bution and should be random. An overall check of 
model adequacy is provided by the Ljung-Box Q 
statistic (Kumar&Anand, 2014). [22] The test sta- 
tistic Q is: 


m » 
r; (e) s 
Qm = n(n+2) ) aR” Xin-r 


and step fourth is Forecasting; once the three 
previous steps of the ARIMA model is over, we can 
obtain forecasted values by estimating the appro- 
priate model free from problems (Padhan, P. C. 
2012). [23]. 
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3.2. GM (1,1) Model 


The meaning of the grey system theory (GST), 
proposed by Deng Julong in 1982 (Ikram et al., 
2019) [24], goes back to the origin of the color 
grey, which is a mixture of white and black. Where 
black indicates a complete lack of information, 
while white indicates complete information,. In 
other words, this theory refers to dealing with the 
lack of information, that is, dealing with small and 
weak data (Wu et al. 2015). [25] In addition, it only 
needs four recent data points and tries to process 
them through the process of grey generation to 
achieve reliable and acceptable accuracy for fu- 
ture prediction (Hu 2017). [26] 

It has become one of the essential methods 
that can be applied to forecast and modelling in 
economic, social, and science fields (Zhou& Ma, 
2020). [27] 

The (GST) is based on the GM (1,1) grey model 
of the first order (i.e. differential equations of the 
first order), where the solution of the differential 
equation is solved with the least square method, 
with one variable (Lin et al. 2011). [28] 

The steps in the construction of the grey fore- 
casting model GM (1,1) are as follows (Tsai et al. 
2017) [29], (Yuan et al 2016) [30]: 

1- Define all of the past data obtained as the 

original series: 


X = {x 1), X(2), ..., XO(M)}, XK) = 0, 


Where X is the non-negative original time se- 
quence (Wang el al 2018) [31], and k = 1, 2,..n. 


2- Perform one primary accumulation gener- 
ating operation sequence (AGO) to add to- 
gether the established original series and 
obtain the following generated series: 


XM = {x (1), XM (2), ..., XY (n)} 
k 
XV(k) = yx (O; b= 1,2)...0. 
i=0 
3- The GM (1, 1) model can be built by estab- 
lishing a first-order differential equation 
with one variable X(t), as follows: 
eet aX) (t)=b, where the parameters 
a is called the development coefficient and is 
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called grey input and b is a constant term (Yang 
el al 2021). [32] 


4- The discrete from of the GM (1, 1) differen- 
tial equation model is expressed as fol- 
lows: 


X(k) + aZ(k) = b; k = 1, 2, 3, ..., 0. 


Where Z® = (Z(2),2(3),....,Z(n)) 
and is said to be mean series of x(k) 


5- In order to obtain the parameters A and B, 
Used the least square method and differ- 
ential and difference equations, that is: 


@= [a,b]’ = (6"B)-*B'Y 


X(2) —Z% (2) 1 

X)(3) —Z%(3) 1 
y=| , Be 

X(n) -ZM(n) 1 


6- we can derive the estimated first-order 
AGO sequenceand the estimated inversed 
AGO sequence , using the least squares 
estimation, and is also called the time re- 
sponse function, and it is as follows (chiou 
et al. 2004)[33]: 


ROH = [xO --| en~at yp 2 = 12,3,....,n 
a a 
7- And then, the equation of time response 


sequence is depicted in the following 
equation: 


RO(k+1) =[XOA)-Alem 42, k= 1,2,3,.....0 

Finally, forecasting equation is generated, and 
the fitted values and predicted values of the se- 
quence X can be calculated. 


KOCK +1) =XV(K41)-XVKE = 


= (x) - =) (1 —e%)e-a 
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Figure 1. Time series plot of electricity consumption in Algeria 
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Source: Author’s calculations using Eviews12 


4. EMPIRICAL RESULTS AND DISCUSSION 


The data used in this study is the annual elec- 
tricity consumption in Algeria. Intended yearly 
data were obtained from https://countryeconomy. 
com/energy-and-environment/electricity-con- 
sumption/algeria, from 1980 to 2021, as presented 
in Figure 1. The dataset consists of a total of 42 
observations. The data set was split from January 
1980 to 2015 for training purposes, whereas the re- 
maining was used for testing purposes from 2016 
to 2021 for both methods (ARIMA and GM (1,1)). 

In Figure 1, the time series plot of electricity 
consumption in Algeria shows it had been increas- 
ing from the starting point to 2006 to 2008. There 
was a slowdown in electricity consumption, but 
again, it increased sharply. This is due to several 
reasons, including population growth on the one 
hand, given that the most significant consumption 
tends to the residential sector, especially in recent 
years, as the state paid attention to infrastructure 
and work to attract foreign investment and pro- 
mote the industrial to obtain sustainable devel- 
opment. (see Table 1) 


Table 1: Results of Descriptive Statistics 
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We note through the above table that the min- 
imum value was 5915, and the maximum value was 
68663; the mean value was positive (26975.38), the 
median value was 20182, and the standard devia- 
tion value was 18877.55. The Electricity consump- 
tion value was positively skewed (0.96284). The 
kurtosis value for electricity consumption is at - 
0.3437. And the count of observations was 42. 


4.1. The ARIMA Model 


To verify the stability, we use the Augmented 
Dicky-Fuller Test (ADF) stability test and the Phil- 
lip-perron Tests (PP). The results of the ADF and 
PP tests are reported in Table 2. 

The variable of electricity consumption is 
non-stationary. The electricity consumption se- 
ries must be differentiated once to be stationary; 
hence, the order of integration d equals 1. 


STANDARD 
MEAN MEDIAN DEVIATION KURTOSIS | SKEWNESS | MINIMUM ) MAXIMUM | COUNT 
26975.38 | 20182 18877.55 -0.3437 0.96284 5915 68663 42 


Source: Author’s calculations using Eviews12 
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Table 2: Augmented Dickey Fuller Unit Root 
and phillip-perron Tests 


Adj. t-Stat Prob.* 

Phillips-Perron test statistic 3.858073 1.0000 
Test critical values: 1% level -4.243644 
5% level -3.544284 
10% level -3.204699 

Adj. t-Stat Prob.* 

Phillips-Perron test statistic -3.833855 0.0267 
Test critical values: 1% level -4.252879 
5% level -3.548490 
10% level -3.207094 

t-Statistic | Prob.* 

Augmented Dickey-Fuller test statistic 3.170168 1.0000 
Test critical values: 1% level -4.243644 
5% level -3.544284 
10% level -3.204699 

Augmented Dickey-Fuller test statistic -3.949801 0.0205 
Test critical values: 1% level ~4.252879 
5% level -3.548490 
10% level -3.207094 


Notes: * indicates statistical significance at the 
1% significance level, ** indicates statistical 
significance at the 5% significance 

Source: Author’s calculations using Eviews12 


Figure 2: Autocorrelation and Partial Autocor- 
relation function. 


Date: 12/31/22 Time: 00:06 

Sample: 1980 2015 

Included observations: 35 
Autocorrelation Partial Correlation AC PAC Q-Stat Prob 

0.547 0.547 11.416 0.001 

0.552 0.361 23.384 0.000 

' 0.343 -0.079 28.142 0.000 
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11 0.060 -0.059 34.088 0.000 
12 0.092 -0.001 34.563 0.001 
13 0.027 -0.045 34.606 0.001 
14 -0.047 -0.189 34.744 0.002 
15 -0.086 -0.010 35.221 0.002 


' 
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' 
' 
' 
' 
' 
' 
! 
| 
' 
' 
' 
' 
' 
' 36.344 0.003 


16 -0.128 -0.011 


Source: Author’s calculations using Eviews12 
Figure 2 below represents the plot of the cor- 
relogram (autocorrelation function, ACF) and 
Autocorrelations (ACF) of the first differenced 
series by lag for lags 1 to 16 of the time series 
of the Electricity consumption in Algeria. 
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Table 3: Model Estimation 


Variable Coefficient Std. Error t-Statistic Prob. 
Cc 1455.241 522.7192 2.783983 0.0089 
AR (1) 0.613647. 0.137115. 4.475423 0.0001 
SIGMASQ 1190862. 280651.8 4.243201 0.0002 
R-squared 0.347306 Meandependentvar 1358.286 
Adjusted R-squared 0.306513 S.D. dependent var 1370.473 
S.E. of regression 1141.274 Akaike info criterion 17.01299 
Sum squared resid 41680169 Schwarz criterion 17.14631 
Log likelihood -294.7274 Hannan-Quinncriter. 17.05901 
F-statistic 8.513797 Durbin-Watson stat 2.461924 
Prob(F-statistic) 0.001085 
Inverted AR Roots 61 


Source: Author’s calculations using Eviews12 


After the experiment, we will select the most 
suitable model for forecasting out of the proposed 
models, with the lowest BIC (Bayesian Information 
Criterion) and AIC (Akaike Information Criterion) 
values. The results indicated that the most appro- 
priate ARIMA model to forecast was ARIMA (1,1,0). 


4.2. Building and Application the GM (1,1) 
Model of Electricity Consumption in Algeria 


Data of electricity consumption in Algeria for 
the past 1980-2015 years were selected as the 
training set, to construct to establish the original 
time series 


xX = {x (4), x2), ..., XO}. 


After that, we calculated the primary accu- 
mulation sequence, to reduce the randomness of 
time series. 

Following, we will GM (1,1) model was estab- 
lished via X1, and to estimate parameters, we use 
the method of least squares. According to data of 
the electricity consumption, we obtain: 

a= — 0.06030557418308459 

b=5563.228281289543 


After obtaining the values of (a) and (0), they 
are introduced into the differential model corre- 
sponding to the GM (1,1) model, and the response 
function could be solved, to obtain the GM (1,1) 
prediction model, as follows: 
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Table 5: Forecasting Values and Performance Evaluation of Actual Values, GM (1,1), ARIMA (1,1,0) for 


Set Testing Sample Data. 


Time Actual value GM (1,1) ARIMA (1,1,0) 
2016 55823 47804.30442 56602.54751 
2017 60283 50775.87094 58057.78893 
2018 62114 539321532 59513.03035 
2019 63695 57284.63334 60968.27178 
2020 66646 60845.50721 62423.51321 

2021 68663 64627.72879 63878.75465 
MAPE 0.679522 0.044765797 


Source: Prepared by researchers 


5§563.228281289543 
—0.06030557418308459 
5563.228281289543 
7 <aassossraisanes) 


RYO(k+1) = [s91s ( )| e(0.0603055741830845%(k~1)) 


5563.228281289543 


Taeelan =9 65° ( Se 
a ae Sige —0,06030557418308459 


)- 92250.65 
Finally, the back difference ratio can be calcu- 
lated, which is: 


S¥ 02 
c= 1/5 Si 
Sy 


Where: 

C is the index of the accuracy of the reaction 
model. 

S? is the variance of the original data: 


1 — 2 
Ss = — Deea(Xp — X) 


S? is the residual al variance: 


1 2 
S? = — Deni (Ex = é) 


Knowing that the smaller the c, the higher the 
accuracy of the model. 
and after testing, 


_2559.0637544640226 


= _ = 0,13720447115478135. 
18651.460356398475 


According to Table 4, which represents accura- 
cy grade reference table. 

C value < 0,35, GM (1,1) The prediction accuracy 
level is: OK or Good. 
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Table 4: Accuracy Grade Reference Table 


Prediction precision grade | Cc 
OK or Good ] < 0.35 
Qualified 0.35 < € < 0.50 
Reluctant 0.50<C < 0.65 
Unqualified C> 0.65 


Source: Zhao, M., Zhao, D., Jiang, Z., Cui, D., 
Li, J., &Shi, X. (2015). [34] The gray prediction 
GM (1, 1) model in traffic forecast applica- 
tion. Mathematical modelling of engineering 
problems, 2(1), 17-22 


Therefore, the GM (1,1) model could be used to 
test the Electricity consumption of Algeria. 

The comparison of the predictive accuracy 
for the above two forecasting models of testing 
sample data is presented in Table 5. The mean ab- 
solute percentage error (MAPE) is utilized, which 
measures the predicting accuracy of the model 
using a Statistical method as defined in Table 1. 
(see Table 5). 

The table above explains that the ARIMA (1,1,0) 
model forecasting is better than GM (1,1) based 
on the mean absolute percentage error (MAPE). 
The MAPE value of testing sample data (2016- 
2021) for GM (1,1) and ARIMA (1,1,0) is 0.679522 and 
0.044765797, respectively. ARIMA (1,1,0) has higher 
forecast accuracy than other forecasting models. 


21 


RIADH KADRI, HACEN KAHOUI, ABDELKADER SAHED 


Figure 2: Curves for the models ARIMA (1,1,0) and GM (1,1). 
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Source: Prepared by researchers 


Figure 2 shows the width of the curves for the 
ARIMA (1,1,0) and GM (1,1) models. The forecasting 
results point out that the values predicted by the 
ARIMA model are very close to the original values 
of Algeria’s electricity consumption from 2016 to 
2020 compared to the GM model. 

Using the ARIMA (1,1,0) model, the annual elec- 
tricity consumption was forecasted until 2030. 

The forecasting values of electricity demand 
are given in Table 6. 


Table 6: Forecasting Value of the ARIMA Model. 


TIME ARIMA (1,1,0) MODEL 
2022 65333.99608 
2023 66789.23751 
2024 68244.47895 
2025 69699.72038 
2026 71154.96182 
2027 72610.20325 
2028 74065.44468 
2029 75520.68612 
2030 76975.92755 


Source: Prepared by researchers 


22 


2020 


Actual value 


GM(1,1) 
» ARIMA(1, 1,0) 


2019 2021 


As revealed by the prediction results forecast- 
ed by the ARIMA (1,1,0) model, based on 2022, elec- 
tricity consumption will continue to increase an- 
nually for the next nine years. 


CONCLUSIONS 


This study aims to provide energy policymak- 
ers and planners with a guideline for choosing the 
most appropriate forecasting techniques to fore- 
cast electricity consumption trends in Algeria ac- 
curately. This paper's main objective is to forecast 
Algeria’s annual electricity consumption for the 
1980-2021 period. Based on GM (1,1) and Autore- 
gressive Integrated Moving Average (ARIMA) to se- 
lect the most appropriate fit. The results show that 
the ARIMA (1,1,0) model exhibits good forecasting 
ability according to the MAPE criteria, Compared 
with Model GM (1,1). 

The Algerian state might use the study’s out- 
come to develop energy policies and measures on 
energy conservation and alternative energy. 

In the future, the following works may focus 
on forecasting electricity consumption using RBF 
Artificial Neural Networks, Hybrid Model GM (1,1)- 
ARCH, and Hybrid ARIMA-GM (1,1). 
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